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PREFACE 


This report presents state-of-the-art irregular wave prediction techniques 
developed by Goda (1975a, 1975b) and suggested by Goda, Takayama, and Suzuki 
(1978). When this report was prepared, little field data were available for 
checking the techniques presented. New data are being gathered at the U.S. 
Army Coastal Engineering Research Center's (CERC) Field Research Facility (FRF) 
at Duck, North Carolina, and results of comparing the new prediction techniques 
with this data are expected in early 1981. 


Much of the design information in the Shore Protection Manual (SPM) assumes 
monochromatic waves; however, input wave conditions for many design analyses 
are often irregular. This report, which supplements SPM Sections 2.3, 3.85, 
and 7.12, presents methods developed for open sections of the coast with contin- 
uously shallowing bottom contours. The major emphasis of the report is on pre- 
dicting nearshore wave height distributions. The work was carried out under 
the offshore breakwater for shore stabilization and wave runup and overtopping 
programs of CERC. 


Goda's methods represent an important increase in the engineering communi- 
ty's ability to predict waves propagating into shallow water. The method is 
based on a number of simplifications and empirical adjustments but appears to 
represent laboratory and limited field data reasonably well. Additional com- 
parisons of predictions with field data are underway. Meanwhile, prediction 
results should be carefully evaluated to confirm that the method is not used 
outside of its range of applicability. 


The report was prepared by William N. Seelig and John P. Ahrens, under the 
general supervision of Dr. R.M. Sorensen, Chief, Coastal Processes and Struc- 
tures Branch. 


Comments on this publication are invited. 


Approved for publication in accordance with Public Law 166, 79th Congress, 
approved 31 July 1945, as supplemented by Public Law 172, 88th Congress, 
approved 7 November 1963. 


TED E. BISHOP 
Colonel, Corps of Engineers 
Commander and Director 
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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to 
metric (SI) units as follows: 


Multiply by To obtain 
inches 25.4 millimeters 
2.54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic -meters 
yards 0.9144 meters 
square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.852 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
millibars 1.0197 x 10°3 kilograms per square centimeter 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.01745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 
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1To obtain Celsius (C) temperature readings from Fahrenheit (F) readings, use 
formula: C = (5/9) (F -32). 


To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 


SYMBOLS AND DEFINITIONS 
stillwater depth 
wave energy moving in direction 0 


wave energy of component ''z"' 


acceleration due to gravity 

deepwater equivalent significant wave height 

deepwater significant wave height 

significant wave height or average of the highest one-third waves 
root-mean-square wave height 

average of the highest l-percent waves 

mean wave height 

constant with energy units 

wave refraction coefficient 

wave refraction coefficient of direction 'Z"' 

deepwater wavelength 

beach slope 

number of direction components used in refraction calculations 
wave setup 

a parameter to quantify the directional spreading of wave energy 
average wave period of the highest one-third waves 

angle of nearshore wave energy vector 

deepwater dominant wave direction 


wave angle 


Ht) 
; i tae ed 


iw tat 


ESTIMATING NEARSHORE CONDITIONS FOR IRREGULAR WAVES 


by 
William W. Seeltg and John P. Ahrens 


I. INTRODUCTION 


Waves are often irregular in height, period, and direction with some or 
most of the waves breaking or near the point of breaking during extreme wave 
conditions. Procedures for predicting design wave conditions for irregular 
waves are not discussed in the Shore Protection Manual (SPM) (U.S. Army, Corps 
of Engineers, Coastal Engineering Research Center, 1977). Much of the infor- 
mation in the SPM is for uniform unidirectional waves. This report supplements 
SPM Sections 2.3, 3.85, and 7.12 with state-of-the-art irregular wave prediction 
techniques developed by Goda (1975a, 1975b) and suggested by Goda, Takayama, 
and Suzuki (1978). Easy-to-use methods for estimating nearshore wave height, 
angle of approach, and resulting water level setup for irregular waves are pre- 
sented. The methods are intended for open sections of the coast with contin- 
uously shallowing depth contours. Wave setup due to narrow band frequency 
spectra and surf beat is considered but other forms of wave-wave and wave- 
current interaction and spectral shape factors have been neglected. Design 
curves and examples of estimating the nearshore significant wave height are 
also available in Seelig (1979). 


The transformation and attenuation of waves propagating from deep water to 
a beach is a problem of considerable difficulty because of a lack of field data 
and a poor theoretical understanding of the complex wave deformation process. 
The methods presented in this report are empirical in nature with the physics 
of the actual problem only partially understood. The results, based on labora- 
tory and limited field data, are considered promising enough to recommend their 
application in selected field calculations. The results of calculations should 
be carefully examined to assure that the basic assumptions of the method have 
not been violated and that reasonable answers result. 


In Section II the directional spreading of a wave energy model suggested by 
Goda, Takayama, and Suzuki (1978) is used to predict the refraction coefficient 
and refracted wave height for the nearshore point of interest. The height is 
then used, in Section III, as input to the surf zone wave height distribution 
model developed by Goda (1975a, 1975b) to estimate the nearshore wave conditions 
and setup. 


The input information necessary for the application of these prediction 
techniques includes the deepwater significant wave height, wave period of peak 
energy density and dominant wave direction, the directional spreading of wave 
energy parameter, S*, beach slope, and water depth at the point of interest. 
In any given design situation, some of these parameters may be more reliable 
than others. Analyses in Section IV show the sensitivity of predicted nearshore 
wave heights to these input parameters. 


Section V applies the techniques presented in Sections II and III to predict 
wave transformation along the 580-meter-long pier at the Coastal Engineering 
Research Center's (CERC) Field Research Facility (FRF) at Duck, North Carolina, 


and compares calculated and measured results. An example is given in Section 
VI that applies the prediction methods to estimate nearshore significant wave 
height. 


II. WAVE REFRACTION - DIRECTIONAL SPREADING OF WAVE ENERGY 


The confused sea state in deep water may be described as the sum of wave 
trains simultaneously moving in various directions. As these wave trains move 
toward the coast, the waves with the largest angles between their crest and the 
bottom contours refract the most, so that nearshore waves appear to be less con- 
fused. The reason for the directional spreading of wave energy in refraction 
calculations for the case of straight parallel bottom contours is discussed in 
this section. 


1. Directional Spreading of Wave Energy. 


Wave direction in deep water is difficult to measure, especially when many 
wave trains of various energy levels are simultaneously moving in different 
directions. However, a few basic measurements will help to quantify this di- 
rectional spreading. Longuet-Higgins, Cartwright, and Smith (1963) suggest the 
following density function for wave energy: 


9.\2 
Eo = (cos) Si (1) 
where 
8 = wave direction angle with respect to the dominant deepwater direction, 
A> (Fig. 1; note that oa 9 is measured from a line perpendicular to 
the shoreline) 
K = a constant used to define the total wave energy 
Eg = the density of wave energy in a given direction, 984 
S* = a parameter that defines the variation of energy level with wave 


direction. 


(Since the frequency dependance of directional spreading of wave energy is of 
secondary importance for diffraction (Goda, Takayama, and Suzuki 1978) and re- 
fraction calculations (Y. Goda, Director, Marine Hydrodynamics Division, Port 

and Harbor Research Institute, Japan, personal communication, 1979) it is neg- 
lected in this report.) Smaller values of S* yield higher amounts of direc- 
tional spreading of wave energy. Goda, Takayama, and Suzuki recommend the values 
of S* in Table 1 for design purposes. 


2. Refraction Calculations. 


Refraction calculations are based on the energy-weighted superposition 
of refraction coefficients obtained from linear theory. If Eg and Kpg are 


Sample values of E,g/K for various combinations ne (eo) tebe | Si 


Se SS 0" 15 30° 60° 90° 
4 1.000 0.934 0.758 0.316 0.063 
12 1.000 0.814 0.435 0.032 0.0003 
37 1.000 0.530 0.077 0.000 0.000 


Energy Distribution With Direction 
Approximated By Breaking The Circle 


{A Into 4,000 Segments, 4 
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Figure 1. Wave refraction including directional spreading 
of wave energy. 
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Table 1. Recommended values of S*. 


S* Wave condition 

4 Wind waves 

2 Swell (short-to-moderate decay distances) 
37 Swell (moderate-to-long decay distances) 


the wave energy and refraction coefficient, respectively, for a wave direction, 
t, then the composite refraction coefficient for waves from several simultaneous 
directions is taken as 


(2) 


where N is the number of wave directions. The computations in Appendix A were 
performed by dividing the 360° circle into 4,000 equal angle segments (N = 4,000) 
and using the spreading function in equation (1). Individual refraction coeffi- 
cients for each segment, Kp;, are computed from the standard methods described 
in Section 2.32 of the SPM or McClenan (1975). The use of equation (2) in the 
calculation is to obtain a better estimate of refraction coefficients for irreg- 
ular waves than would be obtained if a single value of the refraction coefficient 
obtained from linear theory were used. It is recognized that the design curves 
in Appendix A do not account for changes in spectral shape. Note that as S* 
approaches infinity the monochromatic wave refraction solution is obtained 
(agrees with Fig. 2-19 in the SPM). 


3. Wave Refraction Analysis. 


Calculation of refraction coefficients and nearshore wave direction angles 
using the design curves in Appendix A requires the dominant deepwater wave 
direction angle, oa), an estimate of S*, and the representative wave period, 
Tg. For irregular waves the period of peak energy density, which is approxi- 
mately equal to the mean period of the highest one-third waves, should be used 
as the representative wave period. Each figure in Appendix A is for a different 
value of S*. The refraction coefficient, Kp, is determined by finding the 
intersection of the deepwater wave angle, oo, on the abscissa with the value 
of d/(gT3) on the ordinate where d is the nearshore water depth of interest 
and g is the acceleration due to gravity. Kp is estimated by interpolation 
between curves of constant Kp. The nearshore angle of the wave energy vector, 
a (also shown in terms of contours of constant values) is the angle that should 
be used, for example, in determining the component of wave energy flux in the 
longshore direction. 


Table 2 illustrates sample values of Kp and a for various values of 7S* 
and ao where d/ (gTS) = 0.01. ° Note that for small, deepwater wave angles the 
wave refraction coefficient is smaller for lower values of S*; for large deep- 
water wave angles the refraction coefficient is larger for smaller values of S*. 


Table 2. Comparison of refraction coefficients for selected 
S* and a, values at d/(gT3Z) = 0.01. 


Infinity 


(a) 
0° ORI92 OoQ7 0.98 1.00 
ANS 0.82 0.85 0.87 0.88 


0.52 0.44 0.36 


The equivalent deepwater wave height, 4, is determined from 
Ho = Ke Hy (3) 


where Hy, is defined as the deepwater significant wave height. H4 should 
also include diffraction or any other loss coefficients if they are significant. 


4. Example Problem in Wave Refraction. 


GIVEN: The wave period, T, = 10 seconds, a dominant deepwater wave angle, 


> = 40°, and a significant wave height, H, = 2.0 meters. 


FIND: The equivalent deepwater significant wave height, Hj, and the near- 
shore angle, a, for wind waves at a water depth of 1.0 meter. 


SOLUTION: From Table 1 the value of S* = 4 is selected for wind waves. The 
term d/ (gT3) has the value 
1.0 


(9.8 x 10%) 0.001 


From Appendix A, Kp = 0.80 and a = 45 ou so from equation (3) the deepwater 
equivalent wave height is 


HS = Kp Ho = 0.80(2.0) = 1.60 meters 
III. NEARSHORE WAVE BREAKING MODEL 


Goda (1975a, 1975b) developed a nearshore wave height prediction model for 
irregular waves that accounts for wave breaking, nonlinear wave shoaling, ir- 
regular wave setup, and surf beat. A brief description of Goda's model is 
given below. 


1. Characteristics of Goda's Model. 


Goda assumed that the deepwater significant wave height, H,, and the average 
period of the significant waves, T,, are known or can be estimated. The wave 
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heights are assumed to have a Rayleigh distribution in deep water and this dis- 
tribution is used to characterize wave heights as the waves move into shallower 
water until a depth is reached where the waves begin breaking. Goda's approach 
allows the broken waves to reform at a lower height, so that the wave height 
distribution is no longer described by the Rayleigh distribution. As the waves 
move into shallower water the nonlinear method developed by Shuto (1974) is 
used to estimate wave shoaling coefficients. Shuto's method of calculating 
shoaling coefficients usually gives somewhat higher waves than would be pre- 
dicted using the conventional linear shoaling method. Nonlinear shoaling is 
consistent with the observed behavior of waves in shallow water and is conserv- 
ative when compared to linear shoaling. 


The offshore beach profile is assumed to be represented by a straight plane 
surface. Dissipation of wave energy by bottom friction is usually very small 
for typical sand beaches; therefore, bottom dissipation is neglected. 


Wave setup or setdown and surf beat are related to the wave breaking process 
and are accounted for in the model. The radiation stress (i.e., wave-induced 
transport or momentum) of the waves progressing toward the shore causes wave 
setup which can either increase or decrease the local water depth; a decrease 
is often referred to as setdown. Setdown occurs seaward of the breaker zone 
and setup occurs shoreward from the point where a significant number of the 
waves break. Surf beat is the longer period component of water level oscillation 
(periods from 20 seconds to several minutes) due to longer period irregularities 
in wave action. The magnitude of the surf beat is amplified in shallow water. 
Figure 2 shows the conditions used in this model. 


Deep Water 
Effective Water 

Lo 
en d> nom 


Setup, Sw ie 


Setd 
erdown SWL 


Point of Interest 
m, Average bottom slope 
one-half to one wavelength 
seaward of the point of interest 


=| Non- Rayleigh =| Rayleigh 
Ss Ss 
2 3 
a a 
Wave Height Wave Height 
Local Wave Height Distribution Deepwater Wave Height 
Distribution 


Figure 2. Definition sketch. 


2. Design Curves. 


Appendixes B and C present formats of design curves obtained from Goda's 
analytical model (Seelig, 1978). The first format (App. B) includes plots of 
five variables in the nearshore zone: (a) The maximum wave height, Hj , de- 
fined as the mean of the highest 1 percent of the waves; (b) the significant 
wave height, H,, defined as the mean of the highest one-third waves (the 
Significant wave height is approximately equal to four times the root-mean- 
square (rms) surface elevation of a water level record); (c) the rms wave 
height, Hyg; (d) the mean wave height, H; and (e) the wave setup, Sy 
All of these variables are divided by the deepwater significant wave height. 
The five variables are plotted on the ordinate versus the ratio of the local 
Stillwater depth, d, to the deepwater significant wave height, Hee on the 
abscissa. One plot is available for each wave steepness and for each beach 
slope. 


The second format (App. C) gives the ratio of the local significant wave 
height divided by the stillwater depth on the ordinate versus d/(gT?) on the 
abscissa. 


Both formats give the same information about the significant wave height, 
but the first format is often more useful when starting with the deepwater 
parameters H} and Ts. The second format is useful for taking highcast sig- 
nificant wave heights measured at one water depth and using this information to 
estimate the significant height at another shallower depth. 


3. Calculation of Nearshore Wave Heights and Water Level Parameters. 


a. Parameters. Parameters needed to use the curves in Appendixes B and C 
are the equivalent deepwater significant wave height, Hj, the representative 
deepwater wave period, T,, and the offshore bottom slope, m. For an irreg- 
ular profile, take the average profile slope one-half to one wavelength seaward 
of the point of interest. 


b. Procedure for Computation of the Offshore Wave Steepness, H3/L, and 
Depth-to-Height Ratio, d/Hj. If Ht is in feet and ae is in seconds, then 


Ho Hp 


i) Gaeae) we) 


Lj is taken as defined by monochromatic linear theory, so that Lo can be 
considered a reference wavelength. The actual wavelengths during irregular 

wave conditions may be variable and influenced by factors such as spectral 
shape and wave interaction. The advantage of using linear theory to describe 
Lo is that the same wave and setup information can be provided in a more con- 

venient form than if some nonlinear theory were used to characterize wavelength. 
iif HS as) in meters and- 1.) 1s,in)seconds,, then 

0 Hy 

i, | (GL.86 925 > 
O : s 


22) 


The ratio of the water depth at the point of interest to H} is determined by 
using equation (3) 


“c¢. Evaluation of Parameters of Interest. To evaluate the nearshore param- 
eters, use the figure in Appendix B having the offshore slope, m, and offshore 
wave steepness (H3/Lo) closest to the values of interest. Enter the value of 
d/Hj on the ordinate of the graph, and select the curve of the parameter of 
interest. Five dimensionless parameters are presented: 


Curve Description 


1 S/H 

2 Hymg/ Hy, 
3 H/H¢ 

4 H,/H 

5 H/Hg 


Use the abscissa on the right for curve 1 (setup) and the left abscissa for 
curves 2 to 5. Setup has an enlarged scale because it is small compared to 
wave height. 


4. Example Problems to Determine Nearshore Conditions (Neglecting Refraction). 
CeCe St et ee rl Pe CE 2D EXAMPLE PROBLEM 1 * * * * * * * * * * * * * * * & 
GIVEN: The offshore slope is 1/100, HS = 3.0 meters, and T, = 14 seconds. 


FIND: The rms, mean, significant, and maximum wave heights and mean water 
level at a stillwater depth of 6.0 meters. 


SOLUTION: The refracted wave steepness from equation (6) is 


He 3.0 
2. =" = 0.01 
Le A@USol ayia) 


and the depth to deepwater wave height ratio is 


a oe = 2.0 
From Figure B-4.in Appendix B, 
Ratio Value 
Curve Variable value (m) 
1 Siy/HS -0.02 Sy = —006 
2 Hymg/H3 0.87 Hymg = 2.6 
5 H/Hg OVSa yt Hy Wi) ea 
4 Ho/H¢ e2'5 arid Sod 
5 H,/Hg Sy, Hay jen aed 


Mean water level = d + S, = 6.00 - 0.06 = 5.94 meters. 
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ke * k Fk kK RK kK KK kK & * * * EXAMPLE PROBLEM 2% * * * * * ®* * & ¥ ® ® KF RK K 
GIVEN: The average offshore slope is 1/100 and Tg = 8.0 seconds. 


FIND: The variation of significant wave heights at d = 1.50 meters for various 
offshore wave steepnesses assuming no refraction effects. 


SOLUTION: Lo = 1.56 it = 99.8 meters and from Figure C-1 (App. C) for a slope 
of 1 on 100 and d/Lo = 1.5/(1.56(8)7) = 0.015, the following values of H,/d 
and H, are obtained for selected values of H$/lo: 


oh EB Hy Hg at dgyq = 1.5 m 
Lo d 
(m) (m) 
02005 = 025.550.5 0.8 
OAON god ON 10566 1.0 
ONO200 92 Opes O75 1d 
OO ALO) Oni7 il, 2 
OO8t 1/808 WOR85 163 


Depth-limited breaking is important in this example for the larger waves, so 
that as deepwater height increases from 0.5 to 8.0 meters, the nearshore height 
only increases from 0.8 to 1.3 meters in a 1.5-meter water depth. 
CF SF OT 2 PR CC OR ee eer rs AF EXAMPLE PROBLEM 3 * * * * * * * * * * * &® ¥ * * 
GIVEN: A wave gage in 4.0-meter stillwater depth measured a significant wave 
height of 2.0 meters with Tg = 7.6 seconds. The nearshore bottom slope was 
m = 1/50. 


FIND: The significant wave height at a second location with a stillwater depth 
of 1.0 meter. 


SOLUTION: For these conditions at location l, 


Hg wear 
(32), TT Ya 0.50 


d 4.0 
SS O07 
eh (QUSi-x<u7062) 


At location 2, 


(+ eee O ORT Live nons 
ene i (985% 7/562) 


To find the significant wave height at location 2, enter Figure C-2 (m = 1/50) 
for) (d/et-)in =. 0s 007 and find the H3/gT? where (Hg/d); = 0.50. In this case 


ene = 0.0032. Follow the curve for H}/gT? to the point where (d/gT?)o = 
0.0018 and read off the ordinate (Hg/d)» = 0.81. The significant wave height 
at location 2 is 

) x dy 


= @.8 % 1,0 3 Oe 


Note that the ratio of wave height to stillwater depth is large because wave 
setup increases the effective water depth in this case for a steep wave and a 
steep beach. 


IV. SENSITIVITY ANALYSIS 


In any given design situation the quality of the input information may vary. 
An important question tv ask when predicting nearshore wave conditions is: How 
sensitive are the predicted conditions to errors of uncertainty in the input 
parameters? The sensitivity of the predicted nearshore significant wave height 
to input parameters is illustrated below with an example. A reference condition 
is chosen and nearshore significant wave heights are estimated. Each of the 
input parameters is then systematically varied and the results compared to the 
reference condition. 


The reference condition was selected to have the following input parameters: 


H} = deepwater significant wave height (5 meters) 
Ts = wave period (10 seconds) 

= nearshore water depths (0.5 to 30 meters) 
m = beach slope (1/100) 


For the reference condition the predicted nearshore significant wave height is 
slightly less than the deepwater wave height for water depth of 10 to 30 meters 
due to wave setdown (Table 3). The predicted height then decreases in shallow 
water to a value of 0.96 meter in a l-meter water depth (stillwater level). 
Systematic variations are made from the reference condition and from the re- 
sulting predicted nearshore wave heights given in Table 3. Comparisons with 
the reference condition are shown in Figure 3. 


Table 3. Sensitivity analysis for nearshore wave height prediction 
for a selected deepwater condition. 


Stillwater depth (m) 0.5 1.0 3.0 5.0 10.0 20.0 30.0 
Reference d/H} 0.1 0.2 0.6 1.0 2.0 4.0 6.0 


Condition Predicted nearshore 


Reference? 0:7; 0.96 2.1 32 4.9 4.6 4.6 

Variation 
Water 1 m deeper 1.3 1.6 Boll 3.8 5.0 4.7 4.7 
1/20 slope 0.9 eS) 2.6 3.9 Srev2, 4.8 4.8 
Hg =6m ON. 0.9 BE? Soe 5.4 5.6 5.8 
Hg =4mn 0.6 0.9 Ziel Sal 4.0 3.8 ST. 
Ao = AS Ste 4 0.7 0.9 2.0 Sra 4.1 3.9 4.0 
Ty = 658 0.5 0.8 1.9 2.8 beat 4.8 5.0 
Ts = 18 s 0.8 SUS 2.3 Se 6.0 5.2 Sli 


lparameters: Hg = 5.0 meters; Tg = 10.0 seconds; beach slope = 1/100; 
no refraction 


a. Water Depth 1.0m 


deeper than Haan 


bat 
Seceecc 


Yu, 


Figure 3. Sensitivity analysis of nearshore predicted 
significant wave heights. 
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The example in Table 3 and an inspection of Appendixes A, B, and C sug- 
gest the following general rules when predicting nearshore wave height from 
deepwater conditions: 


(a) Water depth becomes an increasingly important factor for 
predictions of wave height in shallower water where (d/H} < 1). 


(b) Deepwater wave height is not important to predictions of 
nearshore wave height where d/H} < 1, but it is more important in 
deeper water where d/H} > 1. 


(c) Substantial errors in wave period produce underestimates 
or overestimates of wave height where d/H} > 1, but are less im- 
portant when d/H} < 1. 


(d) Refraction and directional spreading of wave energy are not 
important variables in shallow areas where d/H} < 1 because refrac- 
tion and breaking effects tend to cancel one another; i.e., ata 
given depth, less refraction is associated with more wave breaking 
and more refraction means less breaking.- Large errors in deepwater 
wave angle will have a pronounced effect where d/H} > 1 or for 
large wave angles. 


(e) In shallow water where d/H} < 1, much of tne information 
about deepwater conditions is lost due to wave breaking and refrac- 
tion effects. 


V. ANALYSIS OF NEARSHORE WAVE HEIGHT CHANGES 
AT THE CERC FIELD RESEARCH FACILITY 


The usefulness of the nearshore wave prediction methods presented in Sec- 
tions II and III is illustrated by comparing observed and predicted nearshore 
wave heights measured at points along the 580-meter-long pier at the CERC 
Field Research Facility (FRF). Note that this is by no means a definitive 
test of the prediction methods for several reasons. First (as discussed in 
Section (IV), some input parameters are more important than others and in the 
case of the FRF data the accuracy of the input information is unknown. Since 
water depth was shown to be especially important in Section IV, profiles were 
taken on either side of the pier and tides were estimated hourly. The pres- 
ence of bars or scour holes under the pier near the gages is not known. Sec- 
ond, the effects of reflection, wave direction, runup, and spray from the 
pier piles on the observed wave records at each of the gages are unknown. 
However, a comparison of observed and predicted wave heights measured at the 
pier is considered useful. 


Since extreme conditions are generally of the most interest to designers, 
waves produced by two storms are examined. The first storm (18 to 23 December 
1977) had a relatively long duration, produced a storm surge of approximately 
0.3 meter, and had maximum significant wave|heights of 3.5 meters at the end 
of pier. A second storm of shorter duration and less wave activity occurred on 
13 September 1978 with maximum significant wave heights of approximately 2 
meters off the end of the pier. 
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Methods developed in this report are illustrated by using the observed 
conditions from one gage to predict wave heights at other locations along the 
pier. For the December 1977 storm, the gage on the end of the pier at station 
580 meters was used to predict significant wave heights at station 190 meters 
(Fig. 4), and observed and predicted waves were compared throughout the storm 
(other gages were not in operation). For the September 1978 storm, a wave 
gage 2,250 meters beyond the end of the pier was used to predict the significant 
wave height at points along the entire pier length. Observed and predicted wave 
conditions were then compared. 


Stotion 
along pier 190 240 325 Wove Gage 30 580 
0 las 
dus. 5 
18 Dec 1977 
10 L141} profile 
0 50 m 
Woter Depth, d (m) Offshore Gage 
ook 2,250 meters 
Survey Date eiotion 1909 240 325 430 580 off end of pier 
lon 
Pree (m) —_——_» 
18 Dec. 1977 1.8 = = = 7.9 
8 Sept. 1978 0.9 4.6 5.1 6.6 8.8 


Figure 4. Profile at the Field Research Facility. 


The beach slope is approximately 1/80 at the FRF, so a diagram of the 
nearshore significant wave height was prepared, similar to the curves in Appen- 
ling (G) (QBs Sys 


The observed tide, wave period of the highest one-third waves, Ts, and 
the significant wave height, H, (defined as four times the standard devia- 
tion of the wave record) are given in Figure 6 for the December 1977 storm. 
Visual observations suggested that the wave angle at the end of the pier was 
small, so a was taken as 0.0°. Since the waves were short-crested (Fig. 7), 
S* = 4 was selected. A prestorm survey indicated that the water depth below 
mean sea level was 7.9 meters at the end of the pier and 1.8 meters at station 
190 where predictions of significant wave height were made using Figures 5 and 
A-3 with the observed conditions at the end of the pier. Observed (solid line) 
and predicted significant wave heights (symbols) show agreement, especially 
during the highest wave conditions (Fig. 6). Tidal variation in water levels 
is associated with wave height changes and this is clearly shown in observed 
and predicted wave heights (not shown in observations at station 190 during 
the end of the storm due to sparse observations). 


For the 13 September 1978 storm, a wave gage 2,250 meters beyond the end 
of the pier was used to record offshore wave height and period. Radar images 
were used to estimate dominant wave direction. Predicted tides and a profile 
survey of 8 September 1978 were used to determine water depth along the pier. 
Figure 8 shows nearshore wave conditions. 
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Figure 5. Nearshore significant wave height diagram for the FRF. 
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Figure 6. Observed and predicted wave conditions at the 
FRF 18 to 22 December 1977. 
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Figure 7a. Waves off the end of the FRF pier (19 December 1977). 
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Figure 7b. Waves along the FRF pier (19 December 1977). 
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Figure 8. Waves at the FRF on 13 September 1978. 
Note waves breaking on the piles. 
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The recommended value of S* = 4 and Figures 5 and A-3 were used to pre- 
dict nearshore wave heights. Predicted and observed wave heights are similar 
for most of the profile (Fig. 9), with the observed wave height in the shal- 
lowest water higher than predicted. 


3.0 
Ho= 1.66 m 
Ts =4.52 5s 
2.0 a,2 50° 
o Predicted Observed 
Hs One 


10.0 


Figure 9. Observed and predicted wave height at the 
FRF, 13 September 1978. 


VI. EXAMPLE PROBLEM (USING NEARSHORE WAVE HEIGHT PREDICTION TECHNIQUES) 


GIVEN: A plane bottom slope of 1 on 100 and the input deepwater wave condi- 
tions in Figure 10. 


FIND: The nearshore significant wave height for water depths between 0.5 and 
6.0 meters for waves coming from each of the three deepwater wave directions. 


SOLUTION: Nearshore conditions are estimated assuming S* = 4 and the deep- 
water wave directions (directions 1, 2, and 3) have absolute wave angles of 
60°, 0°, and 60° with respect to the shoreline. Table 4 summarizes nearshore 
wave predictions for waves from deepwater direction 3; Figure 11 gives wave 
predictions from all directions. 


CONCLUSIONS: Results show that the local significant wave height is primarily 
controlled by depth-limited breaking in 3.0 meters of water or less. In this 
example, wave direction, refraction effects, and wave period are relatively 
unimportant input parameters. The lack of influence of these parameters can 
be seen by examining the predicted nearshore wave heights for directions 1 
and 2. The deepwater wave directions differ by 60°, but the deepwater angle 
has almost no effect on the predicted nearshore wave height in water depths 
of 6.0 meters or less. 


The above example suggests that knowing the total water depth at the point 
of interest is crucial when designing for extreme wave conditions. 


(Xe) 


Direction 2 


Ho= 5.49 m 
Tg? 100s 


Ay: 0° 


Direction 3 Direction 
Ho= 2.41 m Ho=5.70 m 
Ts =6.9 5 T5:10.0 5 
Qo- 60° Qo: 60° 


Shoreline 


Figure 10. Deepwater wave conditions. 


Table 4. Predicted nearshore significant wave 
height (direction 3).} 


dow ade Kr Hy Hi gw Hg 
(m) BT 5 iw Ee oe 
6.0 0.013 0.72 1.74 0.023 3.46 he 72 
4.0 0.0086 0.71 S70) 0.023 C5 S55) Les 
3.0 0.0064 0.70 1.69 0.023 1.76 1.69 
2.0 0.0043 0.70 1.68 0.023 Avealy/ 1.28 
1.0 0.0021 0.69 Wa@7/ 0.922 0.58 0.20 
0.5 0.0011 0.69 1.66 0.022 0.30 0.43 
lH = 2.41 meters, Tg = 6.90 seconds, A, = 60°, 


4. 


Direction 


aN 


0 1.0 2.0 3.0 4.0 5.0 6.0 
dswi (m) 


Figure 11. Predicted nearshore wave height. 
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VII. SUMMARY AND CONCLUSIONS 


Methods developed by Goda (1975a, 1975b) and suggested by Goda, Takayama, 
and Suzuki (1978) for predicting nearshore irregular wave conditions for the 
case of continuously shallowing bottom contours have been presented in forms 
convenient for designers. A sensitivity analysis of the methods shows the 
relative importance of the input parameters on the predicted nearshore sig- 
nificant wave height. Comparison with observed wave height changes at the 
CERC FRF pier indicates that the methods may be suitable for preliminary de- 
sign, although the influence of the pier piles and variable hydrography 
caused by the pier on wave activity has not been determined. Examples are 
presented and design curves are given in Appendixes A, B, and C. 


The design curves in Appendixes B and C were produced by assuming that 
the deepwater wave heights have a Raleigh distribution. In a design situa- 
ation where the deepwater waves are known to be non-Raleigh, which may occur 
with multipeaked spectra, methods described in this report should not be used. 


The analytical model also assumes that the beach has continuously shal- 
lowing bottom contours; the effect of offshore bars on nearshore wave height 
is unknown. As a first approximation for barred coasts the wave height shore- 
ward of the bar should. be taken as equal to the predicted height at the bar 
crest location (Y. Goda, personal communication, 1979). Where the water depth 
becomes less than the depth at the bar crest, the curves can once again be 
used to predict wave heights shoreward of the bar. 


The wave predictions for structural design in very shallow water (approxi- 
mately d/H;, < 0.1) should be used with caution. In very shallow water the 
wave begins to behave like a runup bore with a large forward momentum, and 
for this case the wave height alone may not provide enough information for 
design of a structure. 


When this report was prepared, little wave data were available; therefore, 
only a few observations are presented in Section V. However, large amounts of 
wave data are now becoming available from CERC's Field Research Facility at 
Duck, North Carolina, and CERC is in the process of making a more comprehen- 
sive evaluation of methods described in this report. Preliminary results sug- 
gest that on the average, predictions are conservative for d/H, > 1.5 while 
the predicted significant wave height may be underestimated for d/H, < 1.5. 
Results of CERC's analysis will be published in a later report. 
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APPENDIX A 


WAVE REFRACTION DIAGRAMS 
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Figure A-1. Wave refraction for S* = 4. 
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Figure A-2. Wave refraction for S* = 12. 
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Figure Aq3. 9 Wave) retraction: fom (S* ]= 37. 
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APPENDIX B 


NEARSHORE WAVE CONDITIONS AND SETUP DIAGRAMS 
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Figure C-1. Nearshore significant wave height diagram for 
a 1/100 beach slope. 
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Figure C-2. Nearsnore significant wave height diagram for 
a 1/50 beach slope. 
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Slope m=0.05 (1 on 20) 
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Figure C-3. Nearshore significant wave height diagram for 
a 1/20 beach slope. 
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Figure C-4. Nearshore significant wave height diagram for 
a 1/10 beach slope. 
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